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Abstract 



We show that the dissipative Abelian sandpile on a graph C can be related to a 
random walk on a graph which consists of C extended with a trapping site. From 
this relation it can be shown, using exact results and a scaling assumption, that the 
dissipative sandpiles' correlation length exponent v always equals l/d w , where d w is 
the fractal dimension of the random walker. This leads to a new understanding of the 
known result that v = 1/2 on any Euclidean lattice. Our result is however more general 
and as an example we also present exact data for finite Sierpinski gaskets which fully 
confirm our predictions. 
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Self organised criticality (SOC) [1], |2| is the phenomenon in which a slowly driven system with 
many interacting degrees of freedom evolves spontaneously into a critical state, characterised 
by long range correlations in space and time (for introductory reviews, see H, Q). This 
phenomenon has by now been recognised (or conjectured to exist) in many (models of) 
natural phenomena such as earthquakes ||, forest fires 0, speciation of life 0,. . . Moreover, 
the possible presence of SOC can be investigated in experimentally controlable phenomena 
such as the Barkhausen effect H, rice piles |J, and so on. 

A question which is however still poorly understood is what precisely are the necessary 
ingredients a system (or a model) must have for it to become self-organised critical. For 
example, one may ask whether or not dissipation of 'energy' (or a similar quantity) destroys 
long range correlations. This problem has been studied extensively in the Olami-Feder- 
Christensen (OFC) model of earthquakes 0. Numerical studies originally seemed to show 
convincingly that even in the presence of a small amount of dissipation the model remains 
critical ||10|| . Later, it was shown exactly that at least in mean field, the OFC model is only 



critical when its energy is conserved |TT| . Most recently it was argued on the basis of a study 



of branching rates, that the same is true on a finite dimensional lattice [12 



The situation is less controversial for sandpile models which form the paradigmatic examples 
of systems showing SOC. In the past decade much progress has been made in the theoretical 
understanding of these type of models. This is especially true for the Abelian model |], |j , 
where, following the original work of Dhar a mathematical formalism was developped 
[10-14] that allows an exact calculation of several properties of the model such as height 
probabilities (14|, the upper critical dimension [|18j and so on. 



The role of conservation (of sand) in sandpile models was first studied numerically by Ghaffari 



et al. [19], who found that any amount of dissipation destroys the presence of SOC in the 
model. Recently, it was proven that indeed on any hypercubic lattice a non- conservative 
Abelian sandpile model is not critical ||20|| . There exist therefore a correlation length £ in the 



system which diverges when the dissipation rate goes to zero. This allows the introduction of 
an exponent v which describes this divergence. Numerically fl9|] it was found that v m 1/2 



in d — 2. The same authors argued on the basis of a renormalisation group calculation that 
v = 1/2 on any Euclidean lattice. This result was recently proven exactly f2l| . 



In this paper we study further the Abelian sandpile model with dissipation. We begin by 
showing that this problem can be related to that of a suitably defined random walker on 
a lattice with a trap. This result is quite general and extends an earlier mapping between 
conservative sandpiles and resistor networks (or equivalently the q — > Potts model) [ [Ufl . 
Indeed, we will show that sending the dissipation to zero is equivalent to taking the long 
time limit of the random walk problem. It therefore comes as no surprise that the (sandpile) 
correlation length exponent v can be related to the exponent l/d w which describes the 
asymptotic behaviour of the random walker. Our result implies that v = 1/2 for any 
Euclidean lattice. We thus recover in this situation the conclusion of [IT91 EIH but add a new 



understanding of it. However our prediction is more general and holds also on, for example, 
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fractal lattices, or for certain types of random sandpiles. As an example, we performed 
calculations on the Sierpinski gasket for which d w is known exactly (d w = log 5/ log 2). Our 
data are consistent with the prediction v = l/d w . 

The Abelian sandpile model on an arbitrary graph C (with N vertices) is defined as follows. 
On each vertex i of the graph, there exists a height variable zi which assumes integer values 
and has the interpretation of energy or number of sand grains at site %. The dynamics of the 
model consists of two steps. First, we choose any site i at random, and add one grain of sand 
to that site, i.e. — > Z{ + 1. When at a given site i, the height of sand exceeds a threshold 
z ic , i.e. for Zi > Zi C , we say that that site becomes unstable and then the grains of sand on 
i are distributed among neighbouring sites. This process, called toppling, is specified by a 
matrix A^- such that 

Zj -> Zj — Aij (1) 

The elements satisfy A# > and A^- < when i ^ j, and the condition Ajj > 
which garantees that no sand is created in the toppling process. We will also limit ourselves 
to cases in which A is symmetric. Through toppling, neighbouring sites can become unstable 
and in this way an avalanche of topplings is generated. A new grain of sand is added only 
when the avalanche is over, i.e. all sites are stable again. Finally, grains of sand can leave the 
system through certain boundary sites. These are necessary for the model to reach a steady 
state asymptotically. The number of topplings in a given avalanche s is a random variable 



whose distribution P(s,N) is now known to have rich, multifractal properties ||22|| . In the 
present paper we will however only be interested in the first moment (s) of this distribution. 
An exact expression for this quantity can be obtained as follows [13|]. One first introduces 



the matrix G, which is the inverse of A. The element Gij can be interpreted as the expected 
number of topplings at site j when a grain of sand has been dropped at site i [fl3f1 . From 
this interpretation one obtains 

( S ) = ^EE G* (2) 

One of the results of this paper will be a scaling expression for (s) for the dissipative sandpile 
model. 

To continue we will reason further with the case of a graph where each site (apart from 
the boundary sites) is connected to a fixed number of neighbours z. We take z ic = z,Wi. 
For example, on the square lattice or on the Sierpinski gasket z — 4. In the case of the 
conservative sandpile model we choose the matrix A as 

( z if i = j 

A£ • = < — 1 i and j are neighbours (3) 
[ otherwise 

(in the rest of this paper a superscript c (resp. d) will refer to the conservative (resp. 
dissipative) case). On the other hand, following Tsuchiya and Katori [2D] , we choose in the 
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case of a dissipative sandpile model 



f 27C if i = j 

Af ■ = < — ( i and j are neighbours (4) 
[ otherwise 

with 7 > 1. In this way, at each toppling C z (l ~ 1) grains of sand disappear. 

With each toppling matrix Ay we can associate a random walk problem. To do this we have 
to extend the graph C with one extra site, denoted as T which, as we will see immediately, 
will get the properties of a trap for the random walker. Let C* = C U T. We then define 
a continuous time random walker on C* using the matrix elements Ay as transition rates. 
More concrete, the rate by which the walker jumps from j to i is given by —Ay for any two 
sites on C Secondly, the walker jumps from a site j G £ to the trap T with rate <pj = J2i Ay. 
Once the walker reaches T, it stays there forever. 

Let P(i, k,t) be the conditional probability that the walker is at site i G C* at time t given 
that he was in k at t — 0. This conditional probability then evolves according to the master 
equation 

P(i,k,t) = -J2DijP(j,k,t) (5) 
j 

where 



{Ay i G C, 
-4>j i = T > 
if j = 

For the diagonal elements we take 

/; ,.; = - D ij = - A ij + °j = A n 

In this way, and because of the conditions we put on the matrix A, D has all the necessary 
properties (see e.g. [p3[| ) of a stochastic matrix. 

To solve the master equation (|), it is common (f23fl) to introduce the Green function Gik(s) 
which is the Laplace transform of P(i, k, t) 

G ik (s) = / P(i,k,t)e- St dt (7) 
Jo 



The Green function obeys the linear equation 



3 eC, i^j 

J e C (6) 
T 



Sik - sG ik (s) = D i: jG jk {s) 



whose formal solution is 



This solution can be written in terms of the eigenvalues X a and associated eigenvectors u a 
of D as 

Gik(s) = ( \ A u a )i(u a ) k (9) 

Notice that because of the structure @, the spectrum of the matrix consists of the spectrum 
of A and one zero eigenvalue, which we will denote as Ao- The eigenvector associated with 
this zero eigenvalue is completely concentrated on the trap, (w )j = 8^t- On the other 
hand, for the eigenvectors associated with the other eigenvalues, we have (m q ,) t = 0,a ^ 0. 
Therefore ([J) becomes 

G ik {s) = - s 5 i>T 5 KT + G lk (s) (10) 

where 

Gik(s) = 7 ,\ A Ua)i(u a )k (H) 

We can therefore relate the matrix elements Gij appearing in (||D to the elements of the 
Green function as 

Gy = G y (a = 0) (12) 



We now have all the necessary ingredients to discuss the dissipative sandpiles whose toppling 
matrix is given in ([|). Comparing (0) with (^) , we immediately see that 

A$ = CA& + *C(7 - (13) 

This simple relation implies that the eigenvectors of A d and A c are the same and that their 
eigenvalues are trivially related by a multiplicative and an additive constant. Therefore we 
get for the inverse of Afj using (0) and the general definition ( |Tl~D 

G% = G%{s = 0) = i<%( s = z( 7 - 1)) (14) 

This is our main result. It shows the relation between the dissipative sandpile model and the 
random walker associated with the conservative sandpile model but at s — 2(7 — 1). Taking 
the conservative limit 7 — » 1 then corresponds precisely to taking the limit t — ► 00 in the 
random walk problem. Since that asymptotic limit is determined by the scaling exponent 
d w of the random walk, it can already be expected that the correlation length exponent v is 
related to d w . 

The precise connection between the two exponents can be obtained as follows. We will 
calculate (s) as given in (g) for the dissipative sandpile defined in (p. Summing over i, 
using ([L4p and (^) we get 

1 r°° 

Y.Gt = 7 P %k,t)e-^dt (15) 

T^r C JO 
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where P^k^t) = J2iec P{h k,t) is the probability that the walker that started at k has not 
yet been trapped at time t. To calculate Po(k, t) it is important to first consider the type of 
random walk that we have to investigate. Because of the relation fll^D, we have to work with 
the random walker associated with the matrix A c of the conservative sandpile. The resulting 
random walker is therefore such that only boundary sites, which in the conservative case are 
the only ones where sand leaves the system, are connected with the trap. Random walks 
of this type, at least on an Euclidean lattice are easy to study. We will come back to that 
below. 



First, we consider however the limit iV — > oo. In that limit a random walker starting on 
a typical site will not be trapped for any finite t, since only the boundary sites at infinity 
are connected with the trap. Therefore, one has Po(k, t) = 1. Then the integral in ( |ToD can 
immediately be performed and since the result does not depend on k we get 

(s> = C4^T) (16) 



This relation was first derived on the d = 2 square lattice in [21J but we now see that it is 
quite general. In fact, it shows that the dissipative sandpile is never critical. 



We now turn back to the case that the number of sites on the graph is finite, in which 
case Pq(/c, t) will be a decreasing function of time. Explicit results for this quantity can be 
obtained on Euclidean lattices with elementary Fourier techniques. In one dimension one 
obtains for example 



where u„ 



Similar results can easily be obtained in higher dimensions. In 
the saling limit, L — > oo, t — > oo, it is easy to see that this probability is of the form 

p c (M) ^ (18) 

On any Euclidean lattice, the fractal dimension d w of the type of random walker that we 
consider here equals 2. On the basis of exact results such as flTB|), and on general physical 
intuition it can be expected that for the random walkers which are connected with the trap 
only through some boundary sites, P^k, t) has in general the following scaling behaviour 

^( M )«^ff(£,JL) (19) 

with H some scaling function. Here df is the (fractal) dimension of the graph. 

After inserting (|19|) in fll5l) , making a suitable change of variables, and also performing the 
resulting sum over k we finally get the following scaling form for (s) 

(s) ~L^P(L( 7 - lfl dw ) (20) 
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Figure 1: Plot of (s)L dw versus L(7 - l) 1 /*" for the non conservative sandpile model on 



the Sierpinski gasket with L 



2,4, 



128, and different values of 7. 



where R is a scaling function. 

From (p0|), we see that for the conservative case 
which is an old result obtained by Dhar ||13|| . 



7 = 1, and on an Euclidean lattice (s) ~ L 2 
For 7 > 1, we conclude from (|20|) that the 



Our result ( |20|) is however 
e.g. fractal lattices. As 



exponent v that the describes the crossover between dissipative and conservative sandpiles 
equals l/d w . 

On an Euclidean lattice, we recover in this way a result first determined with an approximate 
renormalisation technique |19[ and recently obtained exactly |H 
much more general. The relation v = l/d w should also hold on 
an example, we checked this scaling on finite Sierpinski gasket. In particular, we considered 
Sierpinski gaskets of n generations with n < 7. For each of these we calculated the matrix Gfj 
by calculating the inverse of Afj using a computer. This we did for several values of 7. From 
G d we then calculated (s) using Finally, we plotted (s)L~ dw versus L(7-1)V*». These 
numerically exact results are shown in figure 1. For the Sierpinski gasket d w = log 5/ log 2. 
As can be seen the agreement with the scaling prediction is excellent. Together with the 
known result for the square lattice case, these data give very strong support for the conjecture 
that v = l/d w . 

One may now ask how general our prediction for v is, and whether one can imagine situations 
in which it does not hold. Besides the scaling assumption (|19D, the crucial step in the 
derivation is the relation (O) between the conservative and dissipative toppling matrices. 
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In a more general case one may consider a toppling matrix of the form 




if % = j 

A? • = — C i and j are neighbours (21) 
otherwise 

where ji is site dependent (the same would hold for a graph in which the coordination 
number is site dependent). Let 7 m = min^. For this case fll4l) gets replaced by 

G% = = z( Jm - 1)) 

where G r is the Green function of a random walker where from site % one enters the trap with 
rate C z (li ~ 7m)- It then depends on the particular values of % what is the precise value of 
d w or even whether d w can still be defined in a meaningful way. It is also not clear in which 
cases a scaling assumption such as (|1^) remains valid. But we expect that for many cases 
that are of interest from the sandpile point of view, our result will hold, eventually with an 
exponent v that depends on the particular distribution of 7^ values. 



As an example, take a matrix Af- whose diagonal elements are constructed in the following 
random way. With probability p we take A^ = zC7ij while with probability 1 — p we take 
Afj = zC72- In the related random walker, and for 72 > 71, the sites that are connected 
with the trap are are randomly distributed over the lattice with probability p. For such a 



random walk, it has been shown that in d — 1, d w — 3 pi |. Hence, we expect that for this 
kind of dissipative sandpile u — 1/3 (instead of 1/2) in d — 1. We are currently verifying 
this prediction. The results will be published elsewhere. 



In summary, we have shown how to relate a dissipative sandpile model with an associated 
random walker. Using results from the theory of random walks we were then able to show 
that v = l/d w for a large class of non-conservative sandpile models. This result is in 
agreement with the available evidence on Euclidean lattices and on the Sierpinski gasket. 
Since the knowledge on random walks is quite extensive, we suspect that many interesting 
phenomena in dissipative sandpiles can now be obtained using the link with random walks. 
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